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1. Introduction 

The relativistic direct interaction theory arises from the expectation that the dy- 
namics of an interacting particle system can be constructed in a consistent Poincare- 
invariant way without introducing the notion of the field as an independent object 
with its own degrees of freedom [ U ]2j, ]3|, f|, |5| . At present the principal possibility 
of such a theory is evident in the classical and quantum domains. Its application to 
the description of particle systems is most effective when processes of radiation and 
particle creation may be neglected. 

Among various more or less equivalent approaches to the construction of the 
relativistic direct interaction theory, the single-time Lagrangian formalism [ [I], [7J 
proposed by Professor Gaida more than twenty years ago, seems to be the most 
convenient for the consideration of the general problem of relativistic dynamics, as 
well as for the investigation of various approximations. This formalism has been 
extended to an arbitrary form of relativistic dynamics [ § defined geometrically by 
means of space-like foliations of the Minkowski space [ |9|, [7], |IOfl . The conditions 
of the Poincare-invariance were reformulated in an arbitrary form of dynamics and 
a wide class of exact solutions to the equations expressing these conditions were 
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established for the interactions originally described by a Fokker-type action. The 
transition from the classical Lagrangian to the Hamiltonian description allows one to 
consider the relativistic effects in the statistical and quantum mechanical properties 
of the particle systems. 

The purpose of the present paper is to review some relatively recent generaliza- 
tions and specifications of this development. The transition from a non-relativistic 
interacting particle system to its relativistic counterpart, which on a more formal 
level can be considered as the replacement of the Galilei group by the Poincare 
group as a symmetry group of the system, leads to profound changes in the struc- 
ture of the theory. Within the Lagrangian formalism such a change manifests itself 
in the necessity of using the interaction Lagrangians depending on derivatives of an 
infinitely high order: in the general case the exact relativistic Lagrangian must be 
defined on the infinite order jet space [ This fact is the Lagrangian counterpart 
of the famous no-interaction theorem in the Hamiltonian relativistic mechanics [ [TT| 



and has with the latter a common cause lying in the very structure of the Poincare 
group. It also reflects the time non-locality inherent to relativistic interactions. All 
the aforementioned exact solutions of Poincare-invariance conditions corresponding 
to time-symmetric Fokker-type actions have such kind of non-locality in any form of 
relativistic mechanics [[12], [13|, fLOfl . Although there are elaborated several methods of 
dealing with such systems (expansions in various parameters [ [I], [HI |15|] , transition 
to the center-of-mass variables [ [H]]), it is evident that such a drastic change in 
the structure of mechanical description leads to serious difficulties in the physical 
interpretation of the formalism, as well as in proving its mathematical consistency. 

But there are important exceptions from the general rule. If the form of dynam- 
ics defines a simultaneity relation in the Poincare-invariant way (i.e. the Poincare 
group transforms simultaneous events into simultaneous ones), then the correspond- 
ing invariance conditions of the Lagrangian description allow a large class of exact 
solutions containing derivatives of any finite order (not less than unity). Particularly, 
in such forms of dynamics we can construct in the closed form a variety of nontrivial 
interaction Lagrangians depending on the first order derivatives. 

This fact was first established for an iV-particle system in the two-dimensional 
space-time M 2 within the framework of the front form of dynamics [ [T7[]. Then it 
was extended to the case of a two-particle system in the four-dimensional Minkowski 
space M4 by means of isotropic forms of dynamics with simultaneity between the 
events of particle world lines defined by a light cone [ |I8| . The existence of such 
"standard" relativistic Lagrangians brings the problem of describing such kind of 
systems within the scope of the usual analytical (and, probably, quantum) mechan- 
ics. It allows the formulation of various exact models of relativistic direct interactions 
which admit more or less explicit investigations. Such models are the main subject 
of this paper. 

It is organized as follows. In section 1 we begin with introducing the notion of 
the form of relativistic dynamics within the framework of the Lagrangian formalism. 
The general features of the relativistic Lagrangian description in a two-dimensional 
version of the front form of dynamics and in isotropic forms of dynamics are discussed 
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in sections 3 and 4, respectively. The Fokker-type action integrals which correspond 
to time-asymmetric interactions are considered in section 5. Section 6 is devoted 
to the construction of the Hamiltonian formalism with constraints in the isotropic 
form of dynamics. On this basis in section 7 we investigate in the most explicit 
form the motions of two particles under the influence of a time-asymmetric scalar, 
vector, and other interactions of physical interest. The limiting case of straight line 
motions of such systems is considered in section 8 within the framework of the front 
form of dynamics. Finally, in section 9 we present certain exactly solvable relativistic 
quantum models of interacting particle systems in the two-dimensional space-time. 



2. Geometrical concept of the forms of dynamics 

Let us consider a dynamical system consisting of iV interacting point parti- 
cles. It is convenient to describe the evolution of this system in the 4-dimensional 
Minkowski space M 4 with coordinates x^, \i = 0, 1, 2, 3. We use the metric \\rj^ u \\ = 
diag(l, — 1, — 1, — 1). The motion of the particles is described by the world lines 
7 a : R — > M4, a = 1, ...,N, which can be parametrized by arbitrary parameters r a . 
In the coordinates we have 

la-Ta^X^Ta)- (2.1) 

The velocity of light is taken to be unity. 

Since in the Poincare-invariant theory no particle can move with the velocity 
greater than the velocity of light, the world lines 7 a must be time-like lines, and the 
tangent vectors 

* = w. (2 - 2) 

obey the inequality 

u 2 a = = u a ■ u a > 0. (2.3) 

It is well known that the whole physical information about the motion of the 
system is contained in the world lines 7 a considered as unparametrized paths in the 
Minkowski space. Therefore, freedom in the choice of parameters r a may be used for 
the simplification of the description. Particularly, we can choose common parameter 
t for all the world lines of the iV-particle system. This parameter is defined by a set 
of iV relations of the following general form: 

<f> a (x 1 (t),...,x N (t),u 1 (t),...,u N (t),t) = 0. (2.4) 

The geometrical concept of the forms of relativistic dynamics originated by Dirac 



[ H [L9| can be introduced within the framework of the single-time Lagrangian or 
Hamiltonian descriptions in the following way [ [| |7j, |TU| . Let us consider the foliation 
£ = {E t \t G R} of the Minkowski space M4 by the hypersurfaces Et with the 
equation 

t = a(x), t E R. (2.5) 
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We require that every hypersurface E t = {x G M n+ i\a(x) = t} must intersect the 
world lines j a of all the particles at one and only one point 

^aW =7aP| S *- ( 2 - 6 ) 



This allows us to consider t as an evolution parameter of the system [ [19|, |20 
In the Poincare-invariant theory, when we consider only time-like world lines, the 
hypersurfaces (2.5) must be space-like or isotropic, 

rj^aWa) > 0, (2.7) 

where <9 M = d/dx^. Then we have d^a > 0, and the hypersurface equation (2.5) has 
the unique solution x° = tp(t, x), where x = (x l ), i = 1,2, 3. Therefore, the constraint 
x a (t) e St enables us to determine the zeroth component of x a (t) in terms of t and 
x' l a (t). The parametric equations (2.1) of the world lines of the particles in the given 
form of dynamics are as follows: 

x° = <p{t, X a (t)) = if a , x i = xl(t). (2.8) 

The evolution of the system is determined by 3N functions t i— > x l a (t). They 
may be considered as representatives (in some local chart) for the sections s : 1R — ► 
F,t h (if:, x l a (t)) of the trivial fibre bundle tt : ¥ — > W with S^-dimensional fibre 
space M. = M. 3N [ |2]J . The latter constitutes the configuration space of our system. 

Three Dirac forms of relativistic dynamics correspond to the following hyper- 
surfaces (2.5): x° = t (instant form), x° — x 3 = t ot x° + x 3 = t (front form), and 
rj^x^x" = t 2 (point form). Other examples may be found in [ [|. 

Now we assume that the evolution of the system under consideration is com- 
pletely determined by specifying the action functional 

S = J dtL. (2.9) 

The Lagrangian function L : J°°7r — ► R is defined on the infinite order jet space of 
the fibre bundle it : F — > IR with the standard coordinates ia S ' [ [22|, |23| . The values 



of these coordinates for the section s : 1 1— > (t,x l a (t)) belonging to the corresponding 
equivalence class from J°°7r are Xa S \t) = d s x^(t)/dt s , s = 0, 1, 2, . . .. The variational 
principle 5S = with the action (2.9) gives Euler-Lagrange equations of motion 



s=0 ® Xa 

where D is an operator of the total time derivative 



dL 

— = °> (2-10) 



OO n o 

D = yy x f^)^ + ^. (2.ii) 
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Let us consider an arbitrary r-parametric Lie group Q acting on M 4 by the point 
transformations g : M4 — > M4: 

x» ^ (gxY = x? + \ a C(%) + o(A), (2.12) 

where X a , a = 1, . . . , r are the parameters of the group. The vector fields 

X a = C<9 M (2.13) 

satisfy the commutation relations of the Lie algebra of group Q, 

[*«, X P\ = cZpXy, a, A 7 = 1, • • • , r, (2.14) 

with the structure constants clg. 

The action (2.12) of group on M4 can be easily extended on the world lines 7 a 
by the rule: 

7a >-> #7a = {fi^k e Im7 a }. (2.15) 

But in the given form of dynamics the world lines 7 a are determined by the functions 
t 1 — > x l a (t) or, in other words, by sections s of the bundle 71". Therefore, (2.15) induces 
an action of group G on J°°7r by the Lie-Backlund transformations [ ^2|, |24], ^3|. As 
it was shown in [ 0], the generators of such transformations have the form: 

^ = ££p s cj— w , (2.16) 

a s=0 °^ Xa 

where 

= Caa ~ v a Vaa, (2-17) 

and 

C = C(*,x a ), V = M(i,x a ), < = < (1) - (2.18) 
The Lie-Backlund vector fields (2.16) obey the same commutation relations as (2.14), 

[X a ,X p ] = cZ fl X 7 , (2.19) 

and commute with the total time derivative (2.11) 

[X a ,D] = 0. (2.20) 

For the Poincare group we have the following ten vector fields corresponding to 
the natural action of V(l, 3) on M 4 : 

< = d„ (2.21) 
X%, = x^d v - x v d» (2.22) 

with the commutation relations 



= o, (2.23) 
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[<, *£] = ^ P X« - ( 2 - 24 ) 

[^1/1 ^p CT ] = Vup^na + Viur^vp ~ Vpp^ua - VucrXftp- (2.25) 

Thus, we obtain the following realization of the Poincare algebra in terms of Lie- 
Backlund vector fields (2.16): 

00 pi 

x l = E E D % - <^\^wv ( 2 - 26 ) 

a s=0 C)Xa 

X i = E E ^'[^^ - - K(x afl a av - x av a a ^}—^, (2.27) 

a s=0 

where we must use (2.8) for the elimination of x° a1 and we denote 

a aix = (d„a)(t,x a ). (2.28) 
Making use of the hypersurface equation (2.5) we find: 

^0 = {dipjdty 1 = <p-\ (2.29) 

0"ai = -¥at(d¥a/dx ai ) EE -^Voi- (2-30) 

It is convenient to introduce the vector fields 

W = — -X^f, Pi = Ji = —-j£ijkXf k , JC = Xm, (2-31) 
obeying the following commutation relations: 

[H, Vi) = 0, \V, n V 3 \ = 0, [H, Ji] = 0, [Pi, J k ) = -e m P h (2.32) 

[JiiJk\ = —tiklJh [K-iiJk] = —ZiklK-h [fc-ijfcj] — ZijkJk-, (2.33) 

[H,Ki]=Pi, [P i ,K j )=6 ij H. (2.34) 

Inserting (2.29), (2.30) into (2.26), (2.27), we obtain the realization of the Poincare 
algebra which is convenient for the consideration of the symmetries of a single-time 
three-dimensional Lagrangian description [ |§ : 

00 Pi 

^ = EE^^-W (2-35) 

OXa 



a s=0 



°° P) 

V * = EE 1 ^' + <^-J]- Wy (2.36) 

a s=Q UX a 



a s=0 

oo 



J-i = e M E ^MO* + <^at)]^L , (2-37) 

^ = E E Vl-fA + - W ai)<Pat]TJM- ^ 

a s=0 ^ Xa 
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The symmetry of the Lagrangian description of an interacting particle system 
under group Q means the invariance of the Euler-Lagrange equation (2.10) under 
corresponding Lie-Backlund transformations generated by the vector fields (2.16). 
The sufficient conditions for the symmetry under the Poincare group have the form 

[0,0: 

X a L = Dtt a , a = 1, . . . , 10, (2.39) 
with auxiliary functions Q a , satisfying the consistency relations 

x a n f ,-Xpn a = cZ f) a r . (2.40) 

An important corollary of symmetry conditions (2.39), (2.40) for an arbitrary 
r-parametric Lie group is the existence of r conservation laws 

DG a = 0, a = l,...,r, (2.41) 

for quantities G a which can be explicitly determined in terms of the Lagrangian 
function L and auxiliary functions Q a . This statement, which is well known as the 
Nother theorem, follows immediately from the identity [ |22|, |24| 



X a L = C aa SaiL + ^i, s D s CL, (2-42) 



which holds for an arbitrary Lie-Backlund vector field (2.16). Here 

dL 



«*, = B-^-^pry (2.43) 



n=s 



are the Ostrogradskyj momenta. Making use of the identity (2.42) in symmetry con- 
ditions (2.39), one readily checks that for the solutions of Euler-Lagrange equation 
(2.10) the conservation laws (2.41) hold with 



oo 



G <* = ^DTaa - ««■ (2.44) 

a s=o 

In the general case the Poincare-invariance conditions forbid the existence of 
interaction Lagrangians which are defined on the jet-space J r ir with some finite r (for 
example, with r = 1). This leads to serious difficulties in the physical interpretation 
of the formalism, and, in fact, makes it impossible to obtain a closed form of the 
corresponding Hamiltonian functions. 

In the following we shall consider some exceptions from this rule. The first is 
offered by the front form of dynamics in the two-dimensional Minkowski space. In 
this case there exists a wide class of interaction Lagrangians for an iV-particle system, 
which are defined on the first-order jet-space J x ix [ The second consists in the 
consideration of a more general definition of the form of dynamics, than (2.5). 
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3. Front form of dynamics in M 2 

In the two-dimensional space-time M 2 the front form of dynamics corresponds 
to the foliation of M 2 by isotropic hyperplanes (i.e., lines): 

x° + x = t. (3.1) 

In this form of dynamics for an iV-particle system Poincare-invariance conditions al- 
low the existence of interaction Lagrangians which do not contain derivatives higher 
than the first order. The general form of such a Lagrangian function including only 
pairwise interactions is given by [ |17| : 

L = - ^m a k a + ^^uKi.hi,^ 1 ,^^ 1 ), (3.2) 

a a < b 



where k a = a/1 — 2v a , r ah = x a — Xb, a, b = 1, N, and V a b are arbitrary functions of 
the indicated arguments. As a result of the Poincare invariance of the Lagrangian 
function (|3.2| ), there exist three conserved quantities: energy E, total momentum P, 
and the center-of-inertia integral of motion K . They have the form [ [L7| : 



dL y> dL 

A' a dv a ' ^ dv a 

a=l a=l 
N 

K = ~tP + J2 x ^ a - ( 3 - 3 ) 



0=1 



The existence of the interaction Lagrangians ( |3.2| ) permits one to trace quite 
easily the relations between various formalisms of relativistic dynamics and to find 
out special features of relativistic particle systems. In spite of the fact that the 
Lagrangian function ( |3.2|) does not contain higher derivatives and the transition to 
the Hamiltonian description is a usual Legendre transformation, the investigation 
of exactly solvable models shows some new features which do not occur in the non- 
relativistic mechanics. 

In the classical mechanics, the Lagrangian function is determined on the tangent 



bundle TAi, L : TAi — > R [ 21 . If the configuration space Ai is diffeomorphic to 



M. N , then the tangent bundle is a trivial one: TAi = M> N x M. N . This means that a 
single chart with coordinates (xx, ...,xjy,Vx, ■■■,Vn) covers the whole TAi. 



For the Lagrangian (|3.2| ) the configuration space Ai coincides with the whole M, N 
or at least with the disconnected union of open sets in R . Hence, one can expect 
that there should not be any complications connected with the global structure. 
But the Lagrangian function ( |3.2|) is determined on submanifold Qf of TAi. This 
submanifold is defined by the inequalities 

v a < 1/2, (3.4) 

which reflect the time-like character of particle world lines in M 2 . Submanifold Qf 
does not have the structure of a tangent bundle. 
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Moreover, we do restrict the Lagrangian description to the smaller region than 
TM. for another reason. The Hamilton principle 5S = leads to Euler-Lagrange 
equations if the Hessian is positively defined: 

h = det\\d 2 L/dv a dv b \ \ > 0. (3.5) 

For the Lagrangian function (|3.2|) the Hessian is, in general, a complicated function 
on coordinate differences and velocities: h = h(r ab ,k c ). Therefore, inequality ( |3.5| ) 
defines an open region Q C TM pa M. 2N . This region also does not have the structure 
of a tangent bundle and for a free-particle system coincides with Qf. 

It could be unimportant if the system moves inside the region (3.5) and does not 
reach the boundary 

dQ = {(x a ,v a ) ETM\h = 0,h- 1 = 0}. (3.6) 

In contrast, the difficulty arises when the system reaches the points of the boundary 
region (singular points) at a finite value of the evolution parameter t [ |||. The 
theorem of existence and uniqueness for Euler-Lagrange differential equations breaks 
at singular points and the Lagrangian system is not defined. Therefore, we cannot 
prolong the evolution of the system beyond the critical points within the framework 
of the basic Lagrangian description. 

The way of overcoming this difficulty is offered by the Hamiltonian descrip- 
tion. It is well known that the Legendre transformation is a differentiable mapping 
£ : TA4 — > T*M.. The transition from the Lagrangian ( |3.2|) to the Hamiltonian 
formalism may be performed by the usual Legendre transformation. But this trans- 
formation is a diffeomorphism only in the region Q. It maps the open region Q C M. 2N 
to the open one £Q C T*M. pa M. A . The Hamiltonian description is equivalent to 
the Lagrangian one only in the region £Q [ ^TJ. In a strict sense the motion in the 
Hamiltonian case is well defined on £Q only. In other words, we should consider 
£Q as a whole phase space of the system. 

After the Legendre transformation is performed, the conserved quantities ( ft.3|) 
become canonical generators of the Poincare group V(l, 1): 

N N 

P+ = ^Pa, K = ^x a p a , (3.7) 

a=l a=l 

N m 2 1 

They satisfy the following Poisson bracket relations: 

{P + ,P_} = 0, {K,P ± } = ±P ± . (3.9) 

Here we have introduced more convenient in the front form quantities P± = E ± P. 
The classical total mass squared function M 2 = P+P- has vanishing Poisson brackets 
with all the generators (|3.7|), fl3.8|). 
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If we deal with the Lagrangian region £Q within the Hamiltonian description, we 
shall obtain the same results as in the Lagrangian case. For systems which reach the 
points of dQ, the Lagrangian description leads to disconnected segments of world 
lines [ [26|] . To obtain the whole evolution of such systems we have to determine 
the motion of the system beyond the Lagrangian region. In the following we shall 
demonstrate for certain relativistic models that the Hamiltonian formalism permits 
one to prolong the evolution of the system beyond singular points and obtain smooth 
world lines in M 2 , as well as in the four- dimensional space-time M 4 (see sections 8 
and 7.1, respectively). 



4. Isotropic forms of dynamics 

For a two-particle system in M 4 the class of isotropic forms of dynamics cor- 
responds to the following definition of simultaneity between the events of particle 
world lines [ |T8[ : 

[xi(t) - x 2 (t)] 2 = (4.1) 
with the supplementary condition 

sgn[x°(t)-x°(t)] 



(4.2) 



where e = ±1. Such a description has been developed within the framework of the 
predictive relativistic mechanics in a series of papers by Kiinzle [ 27, 28, 29]. The idea 
of this definition of simultaneity was formulated in the classic Van Dam-Wigner's 
work [ . In the contents of relativistic Lagrangian and Hamiltonian mechanics the 
descriptions based on equation ( f4.1| ) were elaborated in [ |Tj| |3*I| , 



equations ( |4.1|) , ( f4.2|) determine the difference of the zeroth components: 

z?(t)-^(t)=e| Xl (t)-x 2 (t)|. (4.3) 

For the definition of the value of the common evolution parameter t we choose the 
relation 

'xi(t) + x 2 (t)\ _ ^ 



a 



t. 



where a(x) is the same function as in the definition of the geometrical forms of 
dynamics (|2.5| ). Therefore, we have 



x^t) + x° 2 (t) 



(pit 



Xi(t) +x 2 (t) 



(4.5) 



and 



¥>(*>y) + 2 e l r l 



x° 2 = <p(t,y) - -e|r|. 



(4.6) 



Here and henceforth the variables y M = (x± + x 2 )/2 and r M = x^ — x 2 are used. 
If we put tp(t, y) = t as in the instant form of dynamics, we obtain 



x,. 



e r 



1,2; 



a. 



(4.7) 
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These relations have been used in [ |27|, ^8]. When we choose <r(x) as in the front 
form [ip(t, y) = t — y 3 }, we obtain 

x° a = t + y 3 + ^(-lfe\r\. (4.8) 

In the two-dimensional space-time (|4.8|) reduces to the geometrical definition of the 
front form provided e = sgn(x 2 — £1). 

The general structure of the Lagrange function is again determined by the 
Poincare-invariance conditions. Their formulation requires the realization of alge- 
bra p(l, 3) by the Lie-Backlund vector fields (2.16). In paper it was shown, that 
the components of the corresponding fields have the form (2.17), where 



and 



C = CM*)] (4.9) 



(C a d u a)(^±^-)= Va (t,y)- (4-10) 



All the zeroth components here must be excluded by means of relations ( J4.6[ ). Let 
us note the independence of rj a on the particle labels. 

It is a matter of simple calculation to verify that such vector fields satisfy the 
commutation relations (2.19). 

The Poincare-invariance conditions have the form (2.39) where we can put 

n a = —r] a L. (4.ii) 

Such a choice of auxiliary functions fl a enables (|4.11| ) to be expressed in the form: 

X a L + LDr] a = 0, (4.12) 

where the vector fields 

X a = X a + Va D. (4.13) 

generate the point transformation of the extended configuration space F = R x Ai. 

As in the case of the front form of dynamics in M2, equations ( |4.12j ) allow a 
large class of exact solutions depending on the derivatives of any finite order. If we 
suppose that the Lagrangian contains only the first derivatives, i.e. it is defined on 
the space J 1 ^, we obtain 

dL 2 / dL dL \ 

^dt + E ( CL^r + - < D ^W) + LDr]a = °' (4 ' 14) 

a=l \o a / 



where = C a [x a (t)]. 
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The general solution to these equations can be presented in the form [ [18 



where 



L = 0F(<7i,02,w), 

= (xl - x%)u lfl = (a% - x%)u 2f , = e\r\D(p(t, y) - r • y; 



(4.15) 



r -2 = u^u 



[D<p(t, 



y)-i(-i) a en-v 



n = r/r, r = |r|, v = Vi — v 2 ; 

o"a = = r v u au , u av = u au / a/u 2 ; 



uj = TiT-2 



( J D^(t,y)) 2 -v 1 -v 2 -i(n-v) 2 



u lv u 2 , 



and F being an arbitrary (smooth) function on three variables. 

In the front form of dynamics in M 2 we have d = r, T a = (1 — 2v a 
and uj is a function on the invariants a\, a 2 . 



UJ 



1 + Ea 















V02 





Invariance conditions ( (4.12 ) lead to the conservation laws for the 
(2.44). In our case they have the form: 



dL 



Ga = ^(c-<^)^--a* 



0=1 



-1/2 = k -l 



(4.16) 
quantities 

(4.17) 



Taking into account ( |4.11| ) they can be expressed as 

dL 



Ga - CLa ^ VaH, 



where 



a=l 



Let us introduce the Poincare-invariant functions: 



dF 

A a = &la h ( U(T a ~ °~a 

da a 
dF 



dF 

du ' 

dF 



da„ doj 



(4.18) 
(4.19) 

(4.20) 
(4.21) 



They are not independent, 



axiAt-Bt) = a 2 (A 2 -B 2 ) 



(4.22) 
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In terms of these functions we have: 

dL 1 ~ / OF dF\ „ OF 

— = -( r(pi (t, y) - en)F + v ai T a a a [a a — + ^—)- v bt T- a a a — 

- ^ l (t,y)(T 1 u° l A 1 + T 2 u° 2 A 2 ) + -(-lyemiTmlBt - T 2 u° 2 B 2 ). (4.23) 

where 

* = * < 4 - 24 > 

a=l 

The function ( 4 . 1 9|) is easily found to be 

H = <p t (t, y)(-rF + T.u ^ + T 2 u° 2 A 2 ). (4.25) 
Explicitly, the integrals of motion fl4.18|) are given by 

G a = -£(t, y)<p t (t, yy'H + £ (t, y)(-er,F + T.vuA, + T 2 v 2l A 2 ) + 

+ \{<L - CL)[r 2 (en iW ° - v 2l )B 2 - r 1 (en i «S - v^B,]. (4.26) 

Inserting the expressions for functions which correspond to the generators (2.21), 
(2.22) of the Poincare group, we obtain the following formulae for conserved energy 
E, momentum P, angular momentum J and the center-of- inertia integral of motion 
K: 

E = ip~ l H = -rF + T 1 u < lA 1 + T 2 u° 2 A 2 , (4.27) 
P = - e rF + V 1 v 1 A 1 + r 2 v 2 A 2 , (4.28) 

J = y x P + irx (r lV i£i - T 2 ^ 2 B 2 ), (4.29) 

K = yE- cp{t,y)P - ^[T 2 {ru° 2 - erv 2 )B 2 - T^ru* - erv 1 B 1 ]. (4.30) 

We note that the expressions ( [4.27] ) , ( |4.28| ) can be united into a 4- vector of momen- 
tum Pp, as well as equations (|4.29|) , (r4.30|) represent a 4-tensor of angular momentum 

P^ = er^F - ux/iAx - u 2fl A 2 , (4.31) 

■V = \ {Vi> p n - y» p v ~ r^u^Bx - u 2il B 2 ) - r^uxyBx - u 2v B 2 )) . (4.32) 

Here 

E = -P , Ji = t llk J l \ K i = J oil (4.33) 

and the zeroth components of the 4-vectors x a and u a must be excluded with the 
help of relations (4.6). 

The structure of the motion integrals Q4.31fl , Q4.32D agrees with the results of 
Refs. [ |2^, ESI which were derived within the framework of the predictive relativistic 
dynamics. 

Ten integrals of motion can be used to reduce the integration of equations of 
motion to quadratures. But it is more convenient to preform such a reduction by 
means of the techniques of the constrained Hamiltonian mechanics. 
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5. Fokker-type action and single-time Lagrangians 

One of the possible ways to specify the form of the arbitrary functions entering 
the general solution of the Poincare-invariance conditions is the comparison with the 
Fokker-type relativistic mechanics [ H, U, f£|, the oldest attempt to construct the 
relativistic direct interaction theory which has a relation to the field description. It is 
based on the manifestly Poincare-invariant variational principle formulated in terms 
of four-dimensional coordinates and velocities of the particles. Such a variational 
principle was first introduced for the electromagnetic interaction by Schwarzschild, 
Tetrode, and Fokker at the beginning of this century and developed by various au- 
thors (see Refs. [[1], ||, |33| and references therein). Later this description was extended 



to other relativistic interactions. The equations of motion following from such a vari- 
ational principle explicitly satisfy the demand of relativistic invariance and can be 
compared with the corresponding field theory expressions. However, this approach 
is not free of difficulties both on physical and mathematical levels. The cost for a 
manifestly Poincare-invariant four-dimensional description is the necessity to use 
a many-time formalism which complicates the physical interpretation of its results. 
Mathematically, it is hard to motivate the obtaining of the equations of motion from 
the action integrals which are obviously divergent because the integration is carried 
out on the whole length of the world lines of the particles [ |4| . 

Within the framework of Fokker-type mechanics the dynamics of a relativistic 
particle system is specified in a manifestly Poincare- and reparametrization-invariant 
way on the basis of the variational principle 5S = with the action being given by 

S = Sj — S int , (5.1) 

where 

S f = - ^2,m a I dr aV /uf, (5.2) 

a J 

corresponds to a free-particle system and 



S int = ^2^2 dr a dr b A ab (x a , 

a < b J J 



x b ,u a ,u b ). (5.3) 



determines two-particle interactions. Here A ab are some functions depending on the 
four-dimensional particle coordinates x% and on the first derivatives u%. They have 



the form [ |33|, [35[ : 

A ab = ^u 2 a u 2 b U ab (x a , x b , u a , u b ), (5.4) 

where = u%/ \/uf and function U ab (which we shall call the Fokker potential) 
depends on the following set of the two-body Lorentz scalars: 

Qab = {X a ~ %b) 2 , 0- ab = (X a - X b ) ■ U a , U ab = U a ■ U b \ (5.5) 

that is 

U ab = U ab (g ab , a ab , a ba , u ab ). (5.6) 
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In papers [ [12|, [L3| it was shown that many-time Fokker-type action integrals 
can be transformed into single-time actions with non-local Lagrangians depending 
on the three-dimensional coordinates of the particles and on all the derivatives of 
the coordinates with respect to parameter t. Such Lagrangians provide us with a 
useful tool for the consideration of various approximations [ [12], [13|, [7| , as well as for 



the transition to the predictive relativistic mechanics and Hamiltonian formalism 
| [TJ|, nj. It was demonstrated [ 13| that non-local Lagrangians corresponding to 



the manifestly Poincare-invariant action integrals satisfy the Poincare-invariance 
conditions within the framework of the three-dimensional Lagrangian description of 
interacting particle systems [ ||. The conservation laws which follow from such an 
invariance were investigated via the Nother theorem. Moreover, the non-local single- 
time Lagrangians which are found on the basis of the Fokker-type action integrals 
represent a closed form for a wide class of solutions of equations (2.29) expressing 
the requirements of the invariance of the Lagrangian description of particle systems 
under the Poincare group [ [L3|, |K| . 

If U ab happens to have the special form 

U a b = e a e b u ab S(g ab ), (5.7) 

then action Q5.1| ) describes the electromagnetic interaction of charges e a within the 
framework of the Tetrode-Fokker-Wheeler-Feynman electrodynamics. Such an ap- 
proach has been extended to the interactions which are mediated by massive scalar 
and vector fields [ f|, |33], |35| : 



the scalar case U ab = g a gbG sym (g a b), (5.8) 

the vector case U ab = g a gb^abG sym {g a b)- (5.9) 

In the above, g a is a coupling constant of particle a and G sym (x) = G sym (x 2 ) is a 
time-symmetric Green function of the Klein-Gordon equation 

(□ + K 2 )G sym {x) = An5(x), (5.10) 

where □ = ri^d^d" and k is a mass of the field quanta. Explicitly, we have 

G sym (x) = 5(x 2 ) - &(x 2 ) — ^J!(«V^), (5.11) 

2 V x 2 



where Q(x) is the Heaviside step function and J\(x) is the Bessel function of order 1. 

There exists a wider class of physically important Fokker-type integrals which 
permit a field-theoretical interpretation of interaction between particles. It corre- 
sponds to Fokker potentials of the following form: 

Uab = g a gbf(^ab)G(g ab ), (5.12) 

where f(uj) depends upon the tensor structure of the fields mediating the interaction, 
and G(x) is a symmetrical Green function of the relevant wave equation. In the case 
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of massless fields G(x) = 5(x 2 ). Especially, for interactions mediated by the massless 
field with the given helicity A = ±n we have [ 

f(u) = T n (u). (5.13) 

One more example is a model of confinement interaction [ R7| , for which 



Uab = g a gb<JabCrba5(Qab)- (5.14) 

Equivalently, this model can be presented in the form ( 5.12|) with f(u) = to and the 
Green function G(x) replaced by the "phenomenological propagator" Q(x 2 ). 

Generally, the Fokker-type action with a time-symmetric Green function leads 
to non-local in time Lagrangians and integral- or difference-differential equations of 
motion. It makes the analysis of particle motions a complicated task (except for the 
case of circular motion when the solution may be constructed explicitly [ |38], |39|1). 
An interesting possibility to obtain ordinary differential equations of motion is to 
replace G in the right-hand side of equation ( |5.12| ) by the retarded (advanced) Green 
function of d'Alambert equation 



G e {x) = 20(ex°)(5(x 2 ), e = ±1. (5.15) 

This choice in the case of a two-particle system corresponds to the model with the 
following particle interaction: the advanced field of the first particle acts on the 
second particle and the retarded field of the second particle acts on the first parti- 
cle. Such interactions correspond to the exact solutions of the Poincare-invariance 
conditions considered above in the front and isotropic forms of dynamics [ [1^, 18, |32 



In such models a one-to-one correspondence of points of two particle world lines 
appears naturally, namely, of those points which satisfy the light cone condition: 



r 2 = 0, er u > 0, i.e., er u = |r|, (5.16) 

This correspondence allows one to reduce the Fokker-type integral to a manifestly 
covariant single-time action, 

St= fdT(L + \r 2 ), (5.17) 



where the Lagrangian multiplier A is introduced to take into account condition (|5.16 ) 
as a holonomic constraint (the boundary constraint er° > is also meant). 

An action of this kind occurs when the Fokker potential has a more general 
structure: 

U = f(u,a 1 ,a 2 )G e (r), a x = a 12 , o 2 = o 21 . (5.18) 
The relevant Lagrangian function reads: 

'xi ~ 



L=-^m^xl-^^f, (5.19) 



a=l 



where the dot denotes a derivative on parameter r. It creates a sufficiently wide 
class of two-particle time- asymmetric models. Their study would not be successful 
without an appropriate Hamiltonian description. 
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6. Hamiltonian description in the isotropic form of dynamics 

The Lagrangian description in the configuration space allows a natural tran- 



sition to the manifestly covariant Hamiltonian description with constraints [ [40| |34 
The corresponding phase space T*M4 is a 16-dimensional one with the Poisson 
brackets [...,...]. They have a standard form in terms of covariant coordinates x^ 
and conjugated momenta defined in a usual manner: 

p a „ = dL/dxl (6.1) 



Since the Lagrangian (|5.19 ) and the constraint (|5.16 ) are Poincare- invariant 



there exist ten Nother integrals of motion, 

2 2 



a=l a=l 



In the Hamiltonian description these P M and are generators of the canonical 
realization of the Poincare group. 

By virtue of parametric invariance of the action Q5.171 ), the Lagrangian ( |5.19 ) 



is singular. Hence the canonical Hamiltonian vanishes, while the dynamics of the 
system is determined by the dynamical constraint of the following general form: 

0(P 2 , p 2 ±} P-r,p ± -r)=0, (6.3) 

which appears together with the holonomic constraint ( |5.16j) ; here p± fl = p M — r^P ■ 
p/P ■ r; P^ and p^ = ^(pi, — P2,) are canonical momenta conjugated to y^ i r M , 
respectively. Both the constraints are of the first class, and they unambiguously 
determine the particle dynamics in M4 (i.e. the particle world lines). 

Since no secondary constraints occur, the system possesses 12 physically essential 
degrees of freedom. In order to single them out explicitly, two subsidiary gauge fixing 
constraints are needed. They can be given in the general form: 

X (2/,r,P,p ± ,t) = 0, [x,0]^O, dx/dt ^0. (6.4) 

ij;(y,r,P,p) = 0, [*P,r 2 }^0. (6.5) 

These constraints permit one to eliminate redundant time-like variables x° a and the 
corresponding momenta p a o, and then to pass to the three-dimensional Hamiltonian 
description. 

The gauge fixing constraints do not influence the dynamics of the model, but their 
choice determines specific features of the final description, namely, the reduced phase 
space P (as a submanifold of T*!^), the induced Poisson brackets, and a possible 
choice of variables, in terms of which these brackets take the canonical form. An 
explicit form of observables (i.e. the covariant particle positions, the generators of the 
Poincare group etc.), being functions of the canonical variables of space P, depends 
on a choice of the gauge fixing constraints, too. Thus, using the arbitrariness of this 
choice one can make an effective influence on the structure of the final description. 
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A special choice of the constraint (6.4) in the form 

X = x(y,r,Po,T), (6.6) 



allows one to avoid a well-known no-interaction theorem [[ilj, that is, to pass to such 
a three-dimensional Hamiltonian description of time-asymmetric models in which 
the spatial covariant particle positions x l a (a = 1,2; i = 1,2,3) become canonical 
variables. 

The three-dimensional Hamiltonian description in terms of covariant variables 
is desirable in various aspects. For example, it simplifies the introduction of the in- 
teraction with external fields and allows a position representation on the quantum- 
mechanical level. But this description is not convenient for solving a two-body prob- 
lem, because it does not provide a relevant separation of external and internal degrees 
of freedom. 

Another choice of the gauge fixing constraint (6.4), 

x = y° + tr(A Q dA T /dP ) - r = 0, (6.7) 

where 

Qfiu = r p p u - r u p p , (6.8) 

|P| = v 7 ^, and matrix || A(P/|P|)^|| G 50(1,3), A^P" = 6$\P\ describes the 
Lorentz transformation into the centre-of-mass (CM) reference frame, leads to a 
three-dimensional Hamiltonian description within the framework of the Bakamjian- 
Thomas model [ IT, |2, 43]. Within this description ten generators of the Poincare 



group P M , J pu , as well as the covariant particle positions are the functions of 
canonical variables Q, P, p, 7r. The only arbitrary function entering into expres- 
sions for canonical generators is the total mass \P\= M(p, 7r) of the system which 
determines its internal dynamics. For time-asymmetric models this function is de- 
fined by the mass-shell equation [ [lj|, 32[ which can be derived from the dynamical 



constraint via the following substitution of arguments on the l.-h.s. of (6.3): 

OOOO i i / \ 

P — > M , p±^> — 7r , P-r^-eMp, p±-r^—n-p; here p = \p\ . (6.9) 



Due to the Poincare- invariance of the description, it is sufficient to choose the CM 
reference frame in which P = 0, Q = 0. Accordingly, P = M, J 0i = (i = 1, 2, 3), 
and the components Si = \e^ k Jjk form a 3-vector of the total spin of the system 
(internal angular momentum) which is an integral of motion. At this point the 
problem is reduced to a rotation-invariant problem of some effective single particle 
which is integrable in terms of polar coordinates, 



p = pe p , 7r = 7r p e p + Se^/p. (6.10) 

Here S = |S|; the unit vectors e p , e v are orthogonal to S, they form together with 
S a right-oriented triplet and can be decomposed in terms of the Cartesian unit 
vectors i, j: 

e p = i cosy? +j sin <p, = — isin<£> + j cosy?, (6-11) 
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where <p is a polar angle. 

The corresponding quadratures read: 



t-t =Jdp dTx p {p, M, S)/dM, (6.12) 
<p-<Po = -[dp d7i p (p, M, S)/dS, (6.13) 



where t is an evolution parameter fixed by constraint (6.7) in the CM reference 
frame, and the radial momentum function of p, M, S is defined by the 

mass-shell equation written down in terms of these variables, 

(M 2 , -7T 2 , eMp, -TT • p) = (V, -TT 2 - ^, eMp, = 0. (6.14) 

The solution of the problem given in terms of canonical variables enables one 
to obtain particle world lines in the Minkowski space using the following formulae [ 



TJ, |32|: 



x° a = t + -(-rep, (6.15) 
Particularly, vector r = Xi — x 2 = p characterizes the relative motion of particles. 



7. Time-asymmetric models of particle interactions with long- 
range and confining potentials 

The explicit form of cf) (6.3) depends in a complicated manner on the choice of 
the original Fokker potential. Its construction is the main difficulty which occurs in 
the analysis of time-asymmetric models. Let us split function into two parts: 

0/ + &nt = O, (7.1) 

where 

0/ = \P 2 + ml) + (ml -ml)^+ pi (7.2) 

corresponds to a free-particle system, and 4> int is to be found. Hereafter we refer to 
(pint as the Hamiltonian potential. 

Only few cases are known when function (f>i nt can be constructed explicitly. They 
correspond to the three-parametric Fokker potential 

U = U s + U v + U c = (a s + a v u + a c aia 2 )G e (r) , (7.3) 

where arbitrary constants. The first and the second terms on the r.- 

h.s. of equation (7.3) correspond to the scalar and vector field-type interactions 
with the coupling constants a s and a v , respectively, and the third term describes 
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the confinement interaction (when a c > 0). In the non-relativistic limit this model 
leads to the potential = (a s + ot v )/r + a c r, where r is the distance between the 
particles. 

The corresponding Hamiltonian potential has the form: 

2a s m 1 m 2 + a v (P 2 - mj - ml) ( b±b 2 
(pint = 2a c — - a v 



eP ■ r \eP ■ r 

mim 2 + (bi - a v )ml + (b 2 - a-, 
eP ■ r((b x - a v )(b 2 - a v ) - a%) 



2 2 2a s m 1 m 2 + (bi - a v )m 2 2 + (b 2 - a v )m\ 
-{a s -aj — —777 -^77 — v -iT\ \ 1A ) 



where 

6 o = e(iPT + (-)V-0- ( 7 - 5 ) 

It is worth noting that the interactions are combined in terms of the Hamiltonian 
potential in a non-linear manner. 

For other Fokker potentials approximation methods (such as coupling constant 
expansion) should be applied for the Hamiltonization procedure. Especially for the 
time-asymmetric analogue of the Fokker potential (5.12) the Hamiltonian potential 
in the second order approximation in coupling constant a = g\g 2 reads: 



e jj . r eP ■ r 



where 
and 



h{v) = {{f{v)-\f{p)f-{f{v)f (7.7) 

P 2 — m? — ml , 

u= -. (7.8) 

2m\m 2 

We note that particular cases of (5.12) are the Fokker potentials which correspond 
to the particle interaction via massless linear tensor fields of an arbitrary rank (see 
equation (5.13)) and their superpositions. 

Below we consider some most interesting features of time-asymmetric models 
described in this section. 



7.1. Vector and scalar models 

We begin with vector and scalar time-asymmetric interactions. These models 
are based on the Fokker-type integrals (5.1) with the Fokker potentials U = U v 
and U = U s , respectively (see (7.3)). Both scalar and vector models were partly 
considered earlier (the former in the two-dimensional space-time only) [ |44], |45|, [46|. 



26| , |20| , [28] , [47]]. Our results obtained by means of both the Lagrangian and especially 
the Hamiltonian formulations of these models complete the analysis of their classical 
dynamics. 

The vector and scalar time-asymmetric models present two-body problems lying 
near the border line of those problems the solution of which can be presented in 
a closed form. A lot of analyses can be made analytically. Especially, turning and 
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other important for the integration points are solutions of the third and fourth 
order algebraic equations, while the quadratures cannot be expressed even in terms 
of elliptic and other special functions and, thus, they need computer work. For 
simplicity here we limit ourselves to the case of equal particle rest masses itlq. 

In the non-relativistic limit the vector and scalar interactions reduce to the 
Coulomb interaction with the coupling constant a (namely, a v and a s , respectively). 
Thus, it is convenient to present the specific features of vector and scalar models in 
comparison with the non-relativistic Coulomb system. 

The variety of solutions to the equations of motion of a two-particle system 
consists of a 12-parametric family. The Poincare transformations (which form a 10- 
parametric group) change only the motion of system considered as the whole. This 
motion does not reflect specific features of the models. Here we do not distinguish 
solutions which differ from one another by the Poincare transformations. So, non- 
equivalent solutions form a two-parametric family. It can be parametrized by values 
of total mass M (or energy E in the non-relativistic case) and spin (internal angular 
momentum) S, the pair of integrals of motion. Thus, a variety of all the possible 
solutions is reduced to some subset of (M, S')-plane. We note that parameters m 
and | a | become unessential when using m and r =\a\ /m as units of measurement 
for momentum- and position-like variables, respectively. We also introduce dimen- 
sionless integrals of motion /j, = \Mjm§ and a — S/ \a\. For the convenience we will 
speak about various solutions (namely, phase trajectories, particle trajectories and 
world lines) as if each of them is placed at the corresponding point of (/x, <r)-plane. 

First of all we shall consider a vector model. Qualitatively different types of the 
phase trajectory (three top graphs of figure 1) correspond to three different domains 
T>(+), T>(—) and T>(0) of (//, cr)-plane (the bottom graph of figure 1). 

The number and the position of (//, <r)-domains on (//, cr)-plane are roughly 
in accordance with the non-relativistic case, while the phase trajectories are more 
complicated: they consist of few disconnected branches. It does means that there 
exist few solutions of the Hamiltonian equations of motion at the same values of the 
integrals of motion fi, a. 

Only one branch, namely, 7° for the attraction case (a < 0) and 7I for the 
repulsion case (a > 0), is regular, i.e. it is a relativistic analogue of the phase tra- 
jectory of the Coulomb system and coincides with the latter in a weakly relativistic 
domain of (//, <r)-plane (i.e., fj, ~ 1, and a ^> 1 for the attraction case). If fj, > 1 
(T>(+) domain of (//, <r)-plane), both 7" and 7I exist and correspond to unbounded 
particle trajectories which are analogues to the hyperbolas of the Coulomb prob- 
lem. We note that particle trajectories 7° have a loop-like shape at the points of 
(pL, (j)-plane which are close to /i — 1, a — 1 (figure 2a). This effect becomes more 
evident for bounded states 7° (existing in T>(— )) as the appearance of the perihelion 
advance (figure 2b) (states Y- °f the repulsion case disappear in this domain). In 
the ultrarelativistic case fi — * (the left lower corner of T>(—)) the particles "stick" 
together, so that the distance between them becomes far less than the distance to 
the centre-of-mass (figure 2c). On curve T the regular states correspond to a circu- 
lar motion of particles, so that in domain T>(0) (from below curve J 7 ) the regular 



21 



Isotropic forms of dynamics 



A.Duviryak, V.Shpytko, V.Tretyak 





Figure 1. Vector model. Various domains of (p,,a)- 
plane (bottom graph) and the corresponding 
types of the phase trajectory (three top graphs). 
Curve T on the (/i, <r)-plane is determined by the 
parametric equations: 



a 



(i 



2p 2 + 5p + 4' 

., _ p(2p 2 + 5p + 4) 
2(p + 4) 2 : 

p€ [0,oo[. 
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c) 



Figure 2. Vector model, a < (attrac- 
tion). Regular particle trajectories for 
various values of fx, a. 



a) V{+ 

b) V{- 

c) V{- 



H = 1.01, a = 1.0; 
/x = 0.95, a = 0.68; 
/x = 0.05, a = 0.01. 
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motion is forbidden. 

The remaining branches 7°, 7+, 7+, and 7+ °f the phase trajectory do not have 
non-relativistic counterparts. They exist and are qualitatively similar on the whole 
(//, cr)-plane. These branches present a rather strange motion of particles, so that the 
sign of a does not characterize the interaction as attractive or repulsive. Moreover, 
it turns out natural to sew up the three branches 7+, 7+, and 7+ into a unique 
one 7° (this is shown in figure 1 for the phase trajectory in V(0)), so that the 
resulting motion is as follows: the particles move from an infinite distance between 
them to their collision, go through one another and go away to the distance ~ r , 
draw closer to one another, collide again, and go away to an infinite distance (figure 
3). Branches 7+ and 7" and the corresponding world lines are pathological in the 
sense that the velocities of massive particles tend asymptotically up to the light 
speed. Besides, these solutions contain critical points (namely, collision and turning 
points) in which massive particles reach but not exceed the light speed at a finite 
time. Nevertheless, the particle world lines turn out smooth both at these points and 
everywhere. Another specific feature of the pathological states is that the evolution 
of particles is spread over a semiinfinite interval of the coordinate time while the 
evolution parameter covers the whole real axes (figure 3). 

The scalar model is more intricate, especially for an attractive interaction. There 
are more qualitatively different types of the phase trajectory which correspond to a 
larger number of (//, a")-domains and which consist of more branches (figure 4). 

Among them only one branch is regular, i.e. analogous to the Coulomb phase 
trajectory. It exists in the domains V(l±; 1). Bounded states (in ; 1)) present a 
motion of particles with the perihelion retardance (unlike the advance in the vector 
model). They disappear from below curve J 7 , a > l/Vb, on which the particle 
trajectories become circular. 

In contrast to the case of a vector model, the domain of regular states is bounded 
not only from below, but also from the left where a motion is not forbidden. The 
border lines X + and J, [i > a/5/8 indicate no special changes in the particle 
motion except the appearance of critical points (which corresponds to reaching the 
light speed) on the particle world lines. These provisionally regular states exist in 
the domains X>(1±; 2) and X>(1±;3). The effect of the perihelion retardance grows 
for them (figure 5a), especially in the domain X>(1±; 3); here the particles move as 
if they attract one another at a large distance, while at a small distance ~ r each 
particle repulses another one by a very (but not absolutely) hard core. The particles 
bounce back off this core with the light speed, but their world lines are smooth at this 
critical point (figure 5b). Going to curve J, a < I/a/5, the particle trajectories tend 
(as in the regular case) to circular ones, but in a very strange manner: the particles 
rebound more frequently (figure 5c), so that in the limiting circular trajectories 
(which corresponds to J itself) the set of critical points becomes dense everywhere. 

Apart from the regular or provisionally regular states (which present a reason- 
able behaviour of the particles on the whole) and the pathological ones (which are 
roughly similar to those in the vector model), the attraction (i.e. a < 0) scalar 
model possesses some exotic states which correspond to a bounded particle motion 
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Figure 4. Scalar model. Various do- 
mains of the (/x, (i)-plane (bottom 
graph) and the corresponding types 
of the phase trajectory (six top 
graphs). 

The curves J 7 , J, H±, and on the 
(fx, cr)-plane are defined by the equa- 
tions: 

■p. ,,2 = 27a 

' ^ 2((3+<7 2 ) 3 / 2 +a(9- CT 2 ))' 



J: 



i 



2(l-<7 2 )' 



A± : jU = (Vl + a 2 ± a) /2. 



V5 V3 
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1 




Figure 5. Scalar model, a < (attraction). Various types of bounded particle trajectories. 

a) P(l-; 2) : ji = 0.95, a = 0.68; b) D(l-; 3) : p = 0.9, cr = 0.5; 

c) D(l-; 3) near J : fi = 0.75, cr = 0.3; d) D(2) near H- : n = 0.6, cr = 0.16666; 

e) D(4;2) : /i = 0.3, a = 0.15. 

at a relative distance of order r . These states exist in the domains T>(2), X>(3±), and 
X>(4; 1) — X>(4; 2), i.e., far from the weakly relativistic domain, and thus they have 
no non-relativistic analogues. For example, in the domain T>(2) the particles move 
as if each particle repulses another one by the hard exterior of an empty core inside 
(figure 5d); in the domains ©(4; 1) — D(4; 2) the trajectory of one of the particles 
always lies inside the trajectory of another particle (figure 5e). 

The variety of solutions described above is obtained within the Hamiltonian 
formulation of the vector and scalar models. Within the framework of the La- 
grangian formalism only regular solutions can be reconstructed completely. Besides, 
this framework partially recovers provisionally regular solutions, namely, some seg- 
ments of world lines between the critical points. Other solutions disappear within 
the Lagrangian formalism. 

In our consideration the Lagrangian formalism is primary with respect to the 
Hamiltonian one. Thus, one can conclude at first sight that non-Lagrangian solutions 
have no physical meaning. On the other hand, the Hamiltonian formulation of the 
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models is an important link toward their quantization, and non-Lagrangian solutions 
may contribute to the resulting quantum-mechanical picture. 

These complicated questions are discussed in more detail in sections 8 and 9.1 
where we study the classical and quantum mechanics of the vector and scalar models 
in M 2 . 



7.2. Scalar-vector model 

The purely vector and scalar time-asymmetric models are calculatingly cumber- 
some and present a rather intricate particle dynamics. The case of arbitrary super- 
position of the scalar and vector interaction is not expected to be simpler (though 
it is also solvable). It follows from the complicated structure of the Hamiltonian 
potential (see equation (7.4) with a c — 0). 

In a special case of superposition, a v = na s = k«, k = ±1, the second term of 
(pint (7.4) vanishes. This structure of the dynamical constraint simplifies to a great 
extent the dynamics of the model and makes it similar, in the mathematical respect, 
to the dynamics of a non-relativistic system with the Coulomb interaction. In this 
case one can expect the existence of an additional integral of motion, the relativistic 
analogue of the Runge-Lenz vector. 

Actually, it is easy to guess the structure of this integral of motion working 
within the framework of manifestly covariant Hamiltonian mechanics [ |3T| . For this 
purpose it is convenient to simplify the free-particle term <ftf (10) of the dynamical 
constraint whose cumbersome form obscures the following treatment of the model 
and is caused by a descriptional rather than dynamical reason. Let us perform the 
canonical transformation P M , r M , p M ) i — > (z^, P M , r M , q^), 

% = p»- -\p^ p ^ *" = !/" + -\p^ (> - i-pr**) > ( 7 - 9 ) 

(the variables r M and P^ remain unchanged) . In terms of new variables the dynamical 
constraint takes the form: 

1 9 \, 9 9n (777? — m,) 2 9 a(P 2 — (mi — K777 2 ) 2 ) 

^^4 P 2 (m ' + m ' )+ 4P2 — = ^ ( 7 - 10 ) 

where 

q±^ = P^E^/P ■ r, 'B. vv = T v n v -T v q^ q ± -P = 0. (7.11) 

Then, it is easy to examine that the relativistic analogue of the Runge-Lenz 
vector has the following form: 

^n^ + °< p2 -<;; (7.12) 

where LT^ = 5^ — P a P u / P 2 . It is indeed an integral of motion, i.e. 

[i^,0]«O, [R„r 2 ]=0 (7.13) 
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and satisfies the relations: 

[Rfj,j PA = o, [R„, JxA = -v»\R* + %*R\, (7-14) 

where the Dirac symbol ~ denotes a weak equality. 

The relations (7.14)-(7.15) are similar to those obtained for the Runge-Lenz 
vector of a simple relativistic oscillator and Coulomb models in [ 4"8"fl. These relations 
are essentially nonlinear and thus their group theoretical treatment is complicated. 
In the present paper we limit our study to the case of the CM reference frame in 
which the corresponding Poisson bracket relation can be linearized. 

For this purpose we reformulate (as in the previous cases) the present time- 
asymmetrical model into the framework of the Bakamjian-Thomas model. Then the 
Runge-Lenz vector becomes R^ = (0, R), where 

R = 7T x S + g(M)p/p, (7.16) 

S = p x 7r is a spin of the system, and the total mass satisfies the equation 

d{M) - 7T 2 - 2g(M)/p = 0. (7.17) 

Here 

d(M) = ^ (M 2 - ( mi + m 2 ) 2 ) (M 2 - (m a - m 2 ) 2 ) , (7.18) 

9{M) = m ( m2 _ (mi _ Kni2)2 ) • (7 - 19) 

Besides, in the CM reference frame the covariant particle positions are the following 
functions of the canonical variables: 

, " mi — mi \ 77 



1+ a M2 a )P + e pM' a = 1,2; a = 3 - a. (7.20) 

The Poisson bracket relations for the internal angular momentum (spin) of the 
system S and the Runge-Lenz vector R are similar to those in the non-relativistic 
Coulomb problem: 

{S u Sj} = e^Su, {R t , S 3 } = £ij k R k , {R u Rj} = -d{M)e lJ k S k . (7.21) 

Indeed, when d(M) = 0, equations (7.21) are the relations for generators of the Eu- 
clidian group £(3). In the case d(M) ^ the Si and the normalized R\ = Rij \/\d\ 
generate the group SO (A), when d(M) < 0, and the group SO(l, 3), when d(M) > 0. 
Taking into account equation (7.21) we obtain the following cases for the algebra of 
internal symmetries: 

so(4) for | mi — m 2 \ < M < mi + m 2 , 

e(3) for M = \m± — m 2 \ and M = mi + m 2 , 
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so(l, 3) for < M < \m 1 - m 2 \ and M > m 1 + 



m 2 . 



The existence of the Runge-Lenz vector makes it possible to obtain both the 
relative and particle trajectories traced by vectors p and x a , respectively, without 
an integration. At first we note that these trajectories are flat curves placed on the 
plane orthogonal to the spin of the system, i.e. p ■ S = x a -S = 0. Vector R lies on 
the same plane, i.e. R • S = 0. Multiplying equation (7.16) by p one can obtain the 
relation: 

R-p = gp + S 2 , (7.22) 

where S = |S|. Let (p be an angle between R and p, i.e. Rp = Rp cos (p. Then 
equation (7.22) can be reduced to the canonical equation of a conic section 

p/p = ecosy? — sgng (7.23) 

with the following canonical parameter p and eccentricity e: 



_S 2 _ 2MS 2 _ R I S 2 M 2 — (mi + m 2 ) 2 

^ \g\ \a\ \M 2 — (mi — Km 2 ) 2 \ ' \g\ y a 2 M 2 — (m x — Km 2 ) 2 

(7.24) 

Searching for the equations of particle trajectories is a similar but somewhat 
complicated task. Let us define the vectors: 

r a = x a - c a R, (7.25) 

where 

2(-) s 

(7.26) 



a (M + m a ) 2 -m 2 a ' 
Then, one can obtain the relations 

{-)*R.r a = gr a + ^l5 2 , (7.27) 

which are similar to equation (7.22) and hence can be written down as follows: 

Pa/r a = e cos (f a - sgn g, (7.28) 

where f a are angles between (— ) a R and r a . Equations (7.28) describe the particle 
trajectories as being conic sections of the same shape as the relative trajectory, 
i.e. with the same eccentricity e (7.25) but with other canonical parameters p a = 
jfp. The foci of these conic sections are shifted with respect to the centre of mass 
by vectors c a R. On the contrary, the non-relativistic particle trajectories have a 
common focus which is located in the centre of mass. 
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7.3. Models with higher rank tensor interactions. 

As it was pointed out above, among time-asymmetric field-type models only 
those corresponding to the (arbitrary) superposition of scalar and vector interactions 
permit the exact hamiltonization. In the case when the rank of the field n > 2, the 
transition to the Hamiltonian description and the construction of quadratures can 
be done by means of the method of expansion in a coupling constant. 

The structure of the second order Fokker potential (7.6) is common for linear 
field-type interactions of various tensor dimensions. It specifies the sort of interaction 
by the functions f(u) and h(v) which depend on the integral of motion v only. 
Moreover, the nonlinear gravitational interaction can be also described (at least in 
a slow motion approximation) by this potential (7.6) (see [ f|9], |50|1 ) with 



f gr (u) = 2u 2 - 1, (7.29) 
hg r (u) = -2{2v 2 + 1), (7.30) 

and a gr = -TmiW2 where T is the gravitational constant. It is possible to integrate 
a two-body problem considering / and h as arbitrary first and second order functions, 
respectively. 

We note that in the second order approximation the quadratures for the present 
case can be expressed in terms of elementary functions. For bounded states they 
lead to the relative motion trajectory of a very simple form, 

l/p = \ a \ + 6cos((l -8)<p) (6<H), (7.31) 

where a, b, and 5 are functions of the integrals of motion. It describes an ellipse 
which precesses with the perihelion advance 

Ay? = 2tt5 = -na 2 h{l)/S 2 . (7.32) 

In the case of a linear purely tensor interaction of arbitrary rank n the perihelion 
advance Acp can be calculated by means of the formulae (7.7), (5.12)-(5.13), 

Atp = 7r(2n 2 -l)( gi g 2 /S) 2 . (7.33) 

For the gravitational interaction, using (7.30), we obtain 

Atp = 67r(Tm 1 m 2 /S) 2 . (7.34) 

The spatial particle trajectories calculated by means of (7.16) turn out to be more 
intricate than the relative trajectory which is the typical feature of time-asymmetric 
models. Nevertheless, their analysis leads to the same value of the perihelion advance. 

We note that these relations for the perihelion advance fit those obtained within 
the various quasirelativistic approaches to the relativistic direct interactions [ ET 
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7.4. Confinement models 

Our simplest version of a confinement model [ [55] is based on the Fokker poten- 
tial U c (see equation (7.3)), the time-asymmetric counterpart of which is proposed 
in [ |37H . This model could be regarded as a classical relativisation of the primitive 
quarkonium model with the linear non-relativistic potential. Of course, the relativi- 
sation of any non-relativistic system is not unique. There exists in the literature a 
wide variety of relativistic versions of the potential confinement model. The present 
model has a number of features which are expected for the models of this kind but 
which usually are not realized together. 

1. The model is a self-consistent relativistic two-particle model. The quantities 
in terms of which it is built have a clear physical meaning. Solutions of this model 
are free of any critical point and lead to timelike particle world lines. 

2. It is well known that a non-relativistic potential model with the linear potential 
leads to the Regge trajectory with the unsatisfactory asymptote M ~ S 2 ^ 3 . Here we 
do not propose a quantum version of the present model, but we make the estimates 
of the Regge trajectory from what follows. 

Usually the Regge trajectories in the potential models are calculated in the os- 



cillator approximation [ p6fl. Then, the leading Regge trajectory originates from the 
classical mechanics: it coincides with the curve of circular motions on the (M, S)- 
plane. In our case this curve is described by the following equation: 

(735) 

6v 3a c 

(we consider the case of equal particle rest masses mo). In the ultrarelativistic limit 
M — ► oo it leads to the desirable linear asymptote: 

M 2 « QV3a c S. (7.36) 

It is remarkable that this asymptote is achieved only by taking account of relativity. 

3. The present model permits the interpretation of an interaction in terms of 
some classical fields. It follows from the fact that the Fokker potential U c can be 
transformed into an equivalent form, 

U c = -2a c uD e (x), (7.37) 

where function D e (x), 

D e (x) = ie(ex°)0(x 2 ), (7.38) 
is the fundamental solution of the equation 

n 2 D £ (x) = 4tt5(x). (7.39) 

Thus, the interaction of particles can be considered as mediated by the vector field 
obeying some fourth order equation. Gauge invariant nonlinear equations of this 
kind arise when considering the behaviour of a gluon propagator in the infrared 
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Figure 6. Confinement model. Classical Regge trajectories at various rates of the coupling 
constants and the rest mass. 



region [ [57fl • Static solutions of such equations are used in a sort of the bag model 
of confinement [ [58| . 

The simplest version of the relativistic confinement model can be appropriate for 
the description of light mesons for which the confinement interaction dominates. To 
include into consideration also heavy mesons one can modify the present model by 
adding to U c the usual vector potential U v (with the appropriate coupling constant 



a v < 0) [ 55|. In the non-relativistic limit this mixture leads to the well known 



potential = — |a„| /r + cn c r. The resulting model becomes appreciably cum- 
bersome but still remains solvable. Pathological solutions which occur in this model 
can be unambiguously separated from its regular solutions (which are free of critical 
points). As an illustration we present the classical Regge trajectories for various 
rates of the coupling constants and the rest mass (figure 6). We note that all the 
trajectories tend asymptotically to straight lines. Moreover, the vector correction 
does not influence their asymptotic behaviour which is still described by equation 
(7.36). 



8. Vector and scalar models in M2 

The analysis of the vector and scalar time-asymmetric models in the four-dimen- 
sional space-time M4 was carried out in the previous section for the case of equal 
particle masses. The cumbersome form of the expressions and a large set of possible 
motions obscures the physical understanding of the obtained results. In the two- 
dimensional space-time M 2 the analysis of dynamics becomes considerably simpler 
even for different particle masses. The dynamics in M2 seems to correspond to the 
motions with the inner angular momentum (spin) S = 0. But as it turns out, the 
limit S — > is a singular one. Therefore, the consideration of the dynamics of such 
models in the two-dimensional Minkowski space M 2 appears to be interesting. 

The Fokker-type action integral with a time-asymmetric variant of the Fokker 
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potential (|5.12 ) in the front form in M 2 leads to the Lagrangian [ |17 

JV 



^m a k a 



akik 2 f(u) 



r > 0, 



a=l 



where 



10 



ifhh 
2 U 2 h 



(8.2) 



The existence of three integrals of motion, which for the Lagrangian ( |S.1| ) have the 
form 



K 



mi ^2 olB(uj) 
ki k 2 r 
miki + m 2 k 2 , 



-t(P + +P_)/2-£ 
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' X\k 2 x 2 k\ 



13) 
U) 



a 
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•Ea^la 

k a 

7 2 U 2 h 



X\k 2 x 2 k\ 



i )f 



where 



B(u) = 2(-ujf + (uj 2 -l)f) 



permits one to reduce the solutions of Euler-Lagrange equations to quadrature 
But solutions of Euler-Lagrange equations exist only in the region Q C ThA 
which is defined by the inequalities 



m\m 2 



(8.5) 
(8.6) 

m. 



k 3 k 3 



^(m 2 K, 2 + m\K\)A(u)) ^ 



17) 



where 



A(u) = -f + cof + (uj 2 - l)f" 



The investigation of two-particle models with the time-asymmetric field-like in- 
teractions (see [ |2^, P5| ) shows that for some values of the parameters the system 
reaches the boundary of the Lagrangian region dQ = {(x a ,Xb,v a ,Vb) G M 4 1 = 0; 
hj 1 = 0}. An exception is the repulsion case (a > 0) if the total mass of the system 
M > ni\ + m 2 — Tn, where m 1; m 2 are particle rest masses. Then the system does 
not reach the singular points and the world lines are smooth timelike curves in M 2 
[ [46], [45]. The Hamiltonian description allows one to prolong the evolution of the 
system beyond the critical points for other values of the parameters and, as a result, 
to obtain continuous world lines in the following way 



25 



The Legendre transformation £ associated with the Lagrangian ( |3.1|) with f(uj) 



10 



0,1,2, 



has the form: 



Pa 



m„ 
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dv„ 
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CO 



-i 



(8.9) 
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Here a — 1, 2, a = 3 — a. 

In the scalar (£ = 0) and vector (£ = 1) cases it is possible to solve equations 
with respect to velocities and obtain from the expressions for conserved quantities 
([3.3|) the generators of the Lie algebra of the Poincare group V(l,l) in the explicit 



form [ [)9], [25fl . Separation of the external and internal motions is carried out by the 
choice: 

P+=Pi+p 2 , Q = K/P + - {Q,P + } = 1 (8.10) 
as new external canonical variables. As internal variables we choose 

P + m 
where m = m x + m 2 . Then the Hamiltonian equations of motion become 

M 2 

Q = l/2-—, P + = 0, (8.12) 

• 1 dM " i = -^ 9 -^. (8.13) 



* 2P+ <9£ ' * 2P+ dq 

Solving equations (|8.9|) with respect to velocities and substituting the solutions 
into the expression for the Hessian we obtain from (|8.7p inequalities which define 
the image £Q of the Lagrangian region Q under the Legendre transformation (|8.9| ). 
The external canonical variable P + is an integral of motion. The Hessian does not 
depend on the external variable Q. Thus, all the singularities of the Hessian are 
expressed in terms of inner variables and we can transform ( |8.7| ) into an inequality 
which defines the region £Q in the inner phase space: £Q C M 2 . In the scalar case, 
if a > 0, the region £Q of the phase plane corresponds to the region q > restricted 
by the curves yi, y<i (see figure 7) which are defined by the equations: 

yi : -m^ + m!m 2 + m 2 a/q = , 

(8.14) 

y 2 : m 2 (, + mxm 2 + m^/q = . 

If a < 0, then £Q lies between the curves yi, y 2 to the right of their intersection 
point. 

In the vector case, if a < 0, the region £Q corresponds to the region bounded 
by the curves yi, y 2 , q = (see figure 8) which are defined by the equations: 

?/i : m 1 + £ - a/q = , 

(8.15) 

m ■ m 2 - £ - a/q = . 

If a > 0, then the indicated region lies between the curves y%, y 2 to the right of their 
intersection point. The intersection points of the phase trajectories and the curves 
U2 (yi, 2/2) correspond to the case when one of the particles reaches the speed of 
light: Ki = or k 2 = 0. 
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Figure 7. Scalar interaction. Phase trajectories (continuous curves): (m2—m\)/m = 
0.2; M/m = 1.2. a): a > 0, b): a < 0. Dashed curves yi,y2 correspond to the 
singularity of the Hessian. 
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Figure 8. Vector interaction. Phase trajectories (continuous curves): (m2-mi)/m = 
0.2; M/m = 1.2. a): a > 0, b): a < 0. Dashed curves yi,y2 correspond to the 
singularity of the Hessian. 
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To construct smooth world lines in M 2 it is necessary to consider the inner motion 
in more detail. It is determined by the mass-shell equation 



2 [y 2 — l)m 2 m\m 2 q 2 — 2«M 2 m 1 m 2 m^g+(- l)^ +1 M 4 a 2 



where 
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(8.16) 



(8.17) 
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Figure 9. World lines in M 2 for an unbounded motion: (m 2 — mi)/m = 0.2, M/m = 
1.2 ,a < 0. a): scalar interaction (Stephas case [ ^6|); b): vector interaction (Rudd 
and Hill case [f||). 



We assume that equation ( |8.16 ) is true in the whole phase plane M 2 . The motion 
is possible in the region where 



D e = [v 2 - l)m 2 m(m 2 2 q 2 - 2aM 2 m l m 2 mv l q + (-l) m M 4 a 2 



i.lt 



is non-negative. Then, we see that for a bounded motion q belongs to the interval 
92] j where qi, q 2 are real solutions of the quadratic equation = 0: 



qi 



2aM 2 (-iy +1 
[M 2 — {mi — m 2 ) 2 )m^ 



Q2 



2aM 2 



[M 2 — m 2 )m 



J.19) 



In such a manner we get the phase trajectories which lead to smooth world lines in 
M 2 for all the values of the total mass of the system M > and signs of the coupling 
constant a \ p5|. Using phase trajectory equation (|8.16|) and solving equations QS.12|) , 
( |S.13 ) we obtain a parametric equation for world lines in M 2 : 



xl(q) = t(q) - Xi(q) , x° 2 (q) = t(q) - x 2 (q) ; 



120) 
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Figure 10. Scalar interaction. World lines in M2 for an unbounded motion: (m2 
mi)/m = 0.2, M/m = 1.2 , a < 0. 




Figure 11. Vector interaction. World lines in M2 for an unbounded motion: (rri2 
mi)/m = 0.2, M/m = 1.2 , a < 0. 
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xM = K/P + + (m 2 - £(M 2 , q)) q/P + , 

(8.21) 

x 2 (q) = K/P + - (mi + Z{M\ qj) q/P+ ■ 

Figures 7,8 show examples of the phase trajectories for the scalar (figure 7) and 
vector (figure 8) interactions. Figures 9-11 show the corresponding smooth world 
lines. 

Unlike the scalar interaction, there exist particle collisions in the vector case. 
At the collision points (q = 0) the particles mutually change their positions (figure 
8, b) and the phase trajectories break up. The motion along smooth world lines 
corresponds to the jumps along the momentum axis — oo — > oo (oo — ► — oo). 



9. Quantum models in M2 

In this section we consider a number of exactly solvable quantum-mechanical 
models which follow from certain quantization procedures applied to the correspond- 
ing classical counterparts. We construct a quantum description for the investigated 
above classical time-asymmetric scalar and vector models, as well as for the classical 
models for which the Lagrangian description is not known. 



9.1. Vector and scalar interactions 

The classical two-particle system with time-asymmetric scalar and vector inter- 
actions can be quantized in a purely algebraic way [ |59| regarding the Lie algebra 
so (2, 1) as the basic algebraic structure. Let us introduce the following functions of 
canonical variables: 



Jo 




Jl 


-§( 


J-2 





A q 

.2 Q , a 2 A(a 2 - a 2; 



A 



(9.i; 



where A is an arbitrary constant. They span, under the Poisson bracketing, the Lie 
algebra so (2,1) 

{Jo, J x } = J 2 , {Ji,J 2 } = -Jo, {J 2 ,Jo} = Ji- (9.2) 
Then, the mass-shell equation ( B3B ) takes the form: 

J + C £ = 0. (9.3) 

The quantity 

J = aJ + bJi + dJ 2 (9.4) 
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is an element of the Lie algebra of group SO(2, 1) and we use the following notation: 

M 2 M 2 
a=— — hAm 1 m 2 (m 2 -M 2 ), ^ =_ A — ^ m i m 2( m2 ~M 2 ), 

d=(m2—rrii)(m 2 —M 2 ), Ci=2amm\m2V . (9-5) 

It would appear natural that the structure of the linear relation on the Lie algebra 
so (2,1) must be preserved after quantization. Then, replacing functions fl9.1|) with 
the Hermitian operators obeying the commutation relations of the 50 (2,1) Lie algebra 

[J ,Ji\=U 2 , [Ji,M = -iJo, [J2,M=iJi, (9.6) 

we obtain the quantum-mechanical equation: 

(J + CVM = 0. (9.7) 

This equation was considered in [ |59j as the basic one for the quantum-mechanical 
problem. One can obtain in a purely algebraic way on the basis of equation ( |9.7D 
the mass spectrum 

(M±)| =ml + m 2 2 ± 2m im2 (l - (-1) V/n 2 ) ( ~ 1)V2 , (9.8) 

where 

n = (-1 + v / l + 4(-l)^ 2 )/2 + s, s = l,2,... (9.9) 

The branch (M+) 2 , has a correct non-relativistic limit. Expansion to the order 
1/c 2 gives the following correction to the energy spectrum: 



mim2<y 2 a A m,im2 

E f= 



2ms 2 h 2 Amh A s A c 2 



( 1 _« + -M)i_ 4s( _ 1) , 

s = 1,2,... . (9.10) 



In the single-particle limit (m\/m2 — > 0) we obtain 

/ x(-l) £ /2 

E = mi\l-(-l) i a 2 /n 2 ) - m u (9.11) 

which is in agreement with a one-particle problem in the external scalar or vector 
field in the case of states with the zero value of the quantum orbital number. The 
mass spectrum of a vector type agrees with the result obtained by Barut on the 
basis of the infinite component wave equation [ 60fl. 

The existence of an additional algebraic structure of the mass-shell equation per- 
mits one to quantize the classical problem without ambiguities typical of relativistic 
mechanics [ |61|. Furthermore, such a quantization method allows one to avoid diffi- 
culties connected with the choice of certain representation (coordinate, momentum, 
etc.) which is very important for the field-type interactions because of the difficulties 
of the global structure of the Hamiltonian description (see above). 
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9.2. Relativistic Hamiltonian models in M 2 

Considering the field-type models we started from the Lagrangian description. 
But it is also possible to construct a number of exactly solvable models immediately 



within the framework of the Hamiltonian description [ ^2], ^4|. Contrary to the 
models based on the Fokker-type action integral, relativistic Hamiltonian models are 
not connected with the field theory. Nevertheless, they are also of interest for a vari- 
ety of reasons. They can describe phenomenological aspects of the inner structure of 



mesons and baryons [ |65], |66| . Besides, these models can be useful for the verification 
of different approximation methods, and may be considered as an approximation of 
more realistic models. It appears to be significant for the explanation of relativistic 
effects in the well-established non-relativistic oscillator-like quark models of hadrons. 

The standard quantization procedure consists in the transition from a set of 
canonical generators to a set of Hermitian operators which determine the unitary 
representation of the Poincare group. So, in the case of two-dimensional space-time 
we must put in correspondence with the canonical generators of V(l, 1) the Hermi- 
tian operators K, P + , P_ in some Hilbert space which satisfy the following bracket 
relations: 

[P+,P-]=0, [K,P±] = ±iP ± . (9.12) 
This determines the squared total mass operator M 2 = P+P- and the quantum 



problem is reduced to the eigenvalue problem [ ^2], ^3|, |6~| : 

M 2 ^ = M 2 A V>. (9.13) 

From a variety of the known paths for such a transition we choose the Weyl 
quantization rule [ |67|] . It is necessary that typical of the front form inequalities 

Pa > (9.14) 

be satisfied for this quantization method. It will be noted that these conditions are 
destroyed by field-like interactions. The wave functions ip(p) = (p\ip) describing the 
physical (normalized) states in the front form of dynamics constitute the Hilbert 



space Tifj = £ 2 (M^, dfj,^) with the inner product [ |62], |63], |64|: 



(Vi,V)= / oK(p)V?(?#(p), (9.15) 

where 

N , 

d/4(p) = n ^ e (p«) ( 9 - ie ) 

is a Poincare-invariant measure and Q(p a ) is the Heaviside function. According to 



the Weyl rule we get the following operators [ |62|, |63|, |64 

N N 



P+ = J2p*> K = tJ2Pad/d Pai P_=M 2 /P + , (9.17) 



a=l a=l 
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which are Hermitian with respect to the inner product ( |9.15|) . They determine the 
unitary realization of group P(l.l) on the Hilbert space 7i^. Here M is determined 
by 

M 2 = M 2 f + V, (9.18) 
where M 2 is a free-particle part of the square mass operator: 



N 



^ Pa 



(9.19) 



0=1 



Operator V is an integral operator 



(9.20) 



with the kernel 



V(p,p') 



N 



.d=l 



5(p + - p[) r v (1 

J —00 \ 



Pb+Pb. r lc 



X 



x exp 



N 



^r la (j) a - p' a ) 



a=2 



N 

n 

a=2 



dr la 
2vr ' 



(9.21) 



The general properties of the Weyl transformation [ [67]] ensure that in the classical 
limit these operators correspond to the functions ( |3.7|) , (|3.8|) . 

The evolution of the quantum system is described in the front form of dynamics 
by the Schrodinger-type equation 



where \1/ G Hfj and 



2^ 4 



P-) = -(P + + M 2 /P + ). 



(9.22) 



(9.23) 



Putting \I/ = x(t, P+)ip, where ip is a function of some Poincare-invariant inner 
variables, we obtain a stationary eigenvalue problem for the operator M 2 . In such 
a way a number of exactly solvable two-particle systems were considered in Refs. [ 



62| . p3| . It is convenient to introduce for a two-particle system the following Poincare- 
invariant inner momentum variable f BH1 



V = (Pi -P2V2P+ 



(9.24) 



which is linearly related to the variable £ = (7712 — mi)/2 + mi]. Then the interaction 
part of the squared total mass of the system V takes the form: 



V(r Pl ,rp 2 ) = F(p,r)) 



p = rP + 



(9.25) 
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The conditions (|9.14|) lead to inequality \rj\ < 1/2 . The Hilbert space Ti^ decom- 
poses into the tensor product Ti^ — hi nt <8> Ti-ext ? where "inner" and "external" 
spaces are realized, correspondingly, by functions ip(r]) and x(P+) with the inner 
products 

1/2 

(M) = \ J YJ^2^)4>(V), (9-26) 

-1/2 

oo 

(Xi,X) = I ^XI(P + )X(P+) • (9.27) 
o 

Operator M 2 acts nontrivially only on hi nt . It is an integral operator which is de- 
termined by the rule: 



1/2 



(9.28) 



-1/2 



where kernel W(r], rj') has the form: 

oo 

WM) = 7^JdpF (p, 2±^Q e~^-V) . (9.29) 

— oo 

The structure of operator M 2 coincides with the one-dimensional variant of the 
corresponding expression in Ref. [ |65[| , but in the present treatment kernel W(r),r]') 
is directly related to the classical interaction potential V. 

It is convenient to pass from the functions ip{v) with the inner product Q9.26 ) to 
the functions 

rtv) = ; ^ (9.30) 

with the inner product 

1/2 



(¥>!,¥>)= J d W l( V M V ). (9.31) 

-1/2 

The latter differs from the non-relativistic product only by limits of integration. The 
action of M 2 on the function firj) is defined by the equation 

1/2 

/ 9rn 2 9rn 2 \ f 

(M 2 V )( V ) = + T^J ^ + J ^' W MMv') • (9-32) 

-1/2 
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Let us consider two simple examples. 

1. 5-like potential . Let us put F(p,rj) = aS(p), a = const. Then the equation 



for ip(r]) has the form 



M 



2 



mr, 



1/2 + 7] l/2-r] 



1/2 

) <P{V) = £ f *VW) ■ (9-33) 

-1/2 



Putting 



1/2 



dw(rj) = C & 0) 



-1/2 



we get from ( |9.33|) 



2tt V"" 1/2 + 77 1/2-7? 
Substituting ( |9.35| ) into (|9.34j) , we obtain the equation 

1/2 



-1 



2tt 

a 



-1/2 



, 1 , r ? 771? m 9 

dr] M —± — 

' 1/2 + 7] 1/2-7] 



which describes the eigenvalues of M 2 for bound states. 

/ 



1.5 



-1.5 



-4.5 



\^ /2(A) 


LOG W 2 ' 


2.00 3 



(9.34) 



(9.35) 



(9.36) 



Figure 12. 5-like potential. /i(A) is a graph of the r.-h. side of equation (9.38), /2(A) 



is a graph of the r.-h. side of equation (|9.39[) 
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Let us consider the case of equal particle masses (m\ = m% = m/2). Then we 
get from ( |9.36| ) 

M 

2txM 2 m 2 f dx 



a 2M J x 2 + m 2 - M 2 ' 

—M 



(9.37) 



If M < m, putting M = msinA, < A < 7r/2, we come to the following transcen- 
dental equation for A: 

W ^1 "U-AEAIA). (9.38) 



a \ sin2A^ 

The graph of the right-hand side of this equation (figure 12) shows that there exists 
its only solution for — 3Trm 2 < a < 0. This corresponds to attraction. The energy of 
a bound state has a proper non-relativistic limit. 

It is interesting to point out that there also exists a bound state in the case of a 
strong repulsion. If M > m, one can put M = m ch A, A > 0. Then, from ( |9.37D we 
obtain the following equation: 

W f l+ 2A \ A ^ Xshm (9 . 39) 



a \ sh2A^ 

which has the only solution if a > ixm 2 (figure 12). This solution does not have a 
non-relativistic limit. 

2. Oscillator potential. Let us consider an interaction with a quadratical depen- 
dence on coordinates of the following type: 

V = u 2 r 2 p lP2 = u 2 {l/A - V 2 )p 2 , u e E. (9.40) 

Then, equation ( p. 13 ) transforms into an ordinary differential equation of the hy- 
pergeometric type [ [B2|, 



~ - v 2 ^j <p"(v) - 2 w\v) + 

2 2 \ 1 

mr mi \ 

ip{r,) = (9.41) 



4 + 4 Ml '"' 



2 ' uo 2 y~ n 1/2 + 7] 1/2-7] 
with the boundary conditions 



lim u(7])(p(7]) = 0, lim u(r])(p'(r]) = . (9.42) 

7j — >=bl/2 t; — >=bl/2 



Equation (|9.41|) leads to the mass spectrum 

M 2 = [m + u (n + 1/2)] 2 + ^. (9.43) 

Its nontrivial solutions, which are bounded and square-integrable on the interval 
(-1/2, 1/2) have the form: 

Mv) = c nU + v) U-v) pf "W"°>w-o) (27/) (944) 
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,(2m.2/o;o,2mi/cjo) 



(2r]) are Jacobi polynomials [ ^8| and C n are nor 



In equation fl9.44|) P n 
malization constants. 

In the non-relativistic limit fiwo/mc 2 — *> 0, M n — > m + E n /c 2 we obtain well- 
known wave functions in the momentum representation and a non-relativistic energy 
spectrum of the harmonic oscillator: E = hwo(n + 1/2). 

It is also possible to construct within the framework of the two-dimensional 
variant of the front form an exactly solvable quantum-mechanical iV-particle model 
with the oscillator-like interaction 



(9.45) 



a<b 



Function (|9.45 ) gives an iV-particle generalization of the two-particle interaction 
(|9.40P , as well as one of the possible relativistic generalizations of the N-particle 
oscillator potential. For this system by means of the Weyl quantization rule one can 
also reduce the eigenvalue problem ( |9.13| ) to a differential equation. The system with 
interaction ( |9.45|) has N — 2 additional integrals of motion which mutually commute 
and provide the exact integrability of the system in the classical case. They depend 
nontrivially on the products of coordinate and momenta variables [|64j|. Therefore, 
in general, the quantization procedure can destroy commutation relations between 
these quantities and, as a result, the integrability of the quantum problem. The 
Weyl quantization rule transforms classical additional integrals of motion into a set 
of quantum integrals of motion in involution. That permits one to solve exactly the 
eigenvalue problem and to obtain the eigenfunctions and eigenvalues of M 2 (see [ 



iir = 

n 

Interaction function 
in the coordinates 



N 



N-l 



^rria + ujo ^(rib + 1/2) 



a=l 



6=1 



N 



-UJn 



(9.46) 



]9.45| ) may be generalized by adding terms which are linear 

r ab (Pa-Pb)- (9.47) 



V V = V + a 



a<b 



Such a system also has additional integrals of motion and permits exact solutions 
in the quantum case 



64 



Thus, the Weyl quantization rule preserves the commutation relation of Poin- 
care group V(l, 1), as well as additional symmetries which are responsible for the 
integrability of this model [ |64|. As it was shown in [ ^TJ on the example of a 
two-particle oscillator-like model in the two-dimensional variant of the front form, 
the Weyl quantization is not the only quantization rule with this property. The 
application of different quantization rules preserving the commutation relation of 
V(l, 1) may result in different observables as, for instance, a mass spectrum of the 
system [ plf . 
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10. Conclusion 

We have considered the class of isotropic forms of dynamics which admit the 
construction of a variety of exactly solvable relativistic models of interacting par- 
ticle systems. Most of the models originate from the Fokker-type actions with the 
time-asymmetric (retarded or advanced) Green function of the d'Alembert equation. 
These models reflect not only the relativistic kinematics but also certain field-theory 
aspects of the particle interaction. They demonstrate a complexity of the relativistic 
particle dynamics in comparison with its non-relativistic counterpart. The study of 
such a dynamics in detail is possible because of the fact that the considered forms 
of dynamics allow reformulation of the theory in terms of various formalisms and 
approaches, both three-dimensional and manifestly covariant four-dimensional. 

The physical meaning of time-asymmetric interactions is not so clear. Neverthe- 
less, the corresponding models may be regarded as the first step to some approxima- 
tion scheme for finding solutions of more physically acceptable models, for example, 
the Wheeler- Feynman electrodynamics and the related theories. Particularly, in the 
linear approximation in the coupling constant the time-asymmetric, time-symmetric 
and purely retarded (field) approaches yield the same result. On the other hand, ex- 
act solutions of such models provide a better understanding of the special features of 
relativistic interactions and interrelations between various descriptions of relativistic 
interacting particles. 

We wish to thank Yuriy Kluchkovsky, S. N. Sokolov, V. I. Lengel, and J. Llosa 
for the stimulating discussions. The ideas and influence of the late Professor Roman 
Gaida are evident throughout all the reported investigations. 
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